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Abstract 

Temperature dependence of the quark condensate is studied in the Nambu-Jona- 
Lasinio model with meson loops. Substantial differences are found compared to the 
results with quark loops only. 



1 Introduction 

It is generally believed that QCD undergoes chiral restoration at sufficiently 
high temperatures. This is supported by lattice simulations [1], as well as by 
a variety of model calculations. As the temperature grows, the value of the 
quark condensate increases from its negative T = value and approaches zero. 
As shown in Ref. [2], in the exact chiral limit (zero current quark masses) 
chiral symmetry dictates the form of the first two terms (~ T 2 and T 4 ) 
in the low-temperature expansion of the quark condensate. At higher tem- 
peratures we do not have fundamental knowledge of the behavior of (qq), 
however most of model calculations show a phase transition at temperatures 
T ~ 150 — 200 MeV. Lattice calculations also show a dramatic change of (qq) 
at similar temperatures. 

In this letter we study the temperature dependence of the quark condensate 
in the two-flavor Nambu-Jona-Lasinio (NJL) model [3]. There have already 
been several studies [4-6] of chiral restoration in this model. Our investigation 
brings a new important element: it includes meson loops in a self-consistent 
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way. Previous studies have been performed at the quark-loop level only. At- 
tempts have been made to include meson loops in the NJL model, but self- 
consistency was not completely fulfilled [7]. Our approximation is symmetry- 
conserving [8,9], hence it is consistent with all requirements of chiral symme- 
try. The key ingredient is the self-consistency in solving the equation for the 
scalar density with meson loops present. This makes the approach consistent 
with the requirements of chiral symmetry, such as the Goldstone theorem, 
Gell-Mann-Oaks-Renner and Goldberger-Treiman relations, or one-loop chi- 
ral expansions. 

We find important qualitative and quantitative differences in the tempera- 
ture dependence of the quark condensate in our calculation with meson loops 
compared to the case with quark loops only. With quark loops only, at low 
temperatures the condensate remains flat, whereas in our case it changes con- 
siderably. We show that in the exact chiral limit the change agrees with the 
prediction of the chiral perturbation theory [2] . We also find that meson loops 
decrease the temperature of chiral restoration by about 10%. 



2 Definition of the model 

The Lagrangian of the two-flavor NJL model with scalar-isoscalar and pseu- 
doscalar-isovector interactions is 

£ = q(id^ -m)q + ^ ((qq) 2 + (gi 75 rg) 2 ) , (1) 

where q is the quark field, m is the current quark mass, and 1/a 2 is the 
coupling constant. It is convenient to apply the formalism of effective action 
[10] to Lagrangian (1). Details of this procedure are given in Ref. [9]. Meson 
fields are introduced in the usual way (partial bosonization), with <3> = (<E> , 3>) 
related to the sigma and the pion mean field. At the quark-loop level the 
effective action is 

/($) = J d 4 x (±a 2 $ 2 - a 2 m$ + \a 2 m 2 ) - ±Tr \n(D*D) , (2) 

where D is the Dirac operator , D = d T — ia ■ V + (3& + i/^r • <&. We work 
in Euclidean space-time (r, x). In Eq. (2) we have replaced the usual Tr InD 
term with |Tr In(D^D), which is allowed in the absence of anomalies. In fact, 
this replacement is necessary for the introduction of the proper-time regulator 
[11] used in many NJL calculations, and also in this paper. 
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Meson loops bring an additional term to the effective action [9,10] 



r($) = /($) + |Tr ln^ 1 ) 



(3) 



The inverse meson propagator matrix K is defined as K a ^(x, y) = s<s>a(x)M? b ( y ) ■ 
In Eqs. (2,3) Tr denotes the full trace, including functional space, isospin, 
and in addition color and spinor trace for quarks. In the iV c -counting scheme, 
the quark loop term /($) is the leading contribution of order 0(N C ), and 
the meson loop term |Tr \nK is of order 0(1). Thus the one-meson-loop 
contributions give the first correction to the leading- N c results. 

Using standard methods, Green's functions can be obtained from Eq. (3) via 
differentiation with respect to mean meson fields. Of particular importance is 
the one-point function, which gives the expectation value of the sigma field. 
The condition 



ar($) 



l*o(*)=s 



(4) 



<5>o(x)=S 



yields the equation for the vacuum expectation value of $ , which we denote 
by S. Introducing 



K a (S,Q 2 ) = (4N c f(S,Q 2 )(Q 2 + 4S 2 )+a 2 m/S) 1 , 
K 7r (S,Q 2 ) = (AN c f(S,Q 2 )Q 2 + a 2 m/Sy 1 , 



(5) 



and retaining terms up to order 0(N®), Eq. (4) can be written in the form [9] 



a 2 (S -m) -8N c Sg(S) 



+^/rf 4 Q| 2/(5,0) + ^ (f(S, Q 2 ){Q 2 + AS 2 )) 

K 7T (S,Q 2 ))=0. 



K a (S,Q 2 ) 



+ 3 



d 



2f(S,0) + -^f(S,Q 2 )Q 2 



(6) 



Functions g and / in the above expressions are the quark bubble functions. 
Their form is very simple if no cut-offs were present. In this case we would have 
9(S) = I (0 ^ and f(S, Q 2 ) = J ^ ^ (fc+Q j a+ga , and Eq. (6) could be 
interpreted via standard Feynman diagrams (see Fig. 1). In the presence of a 
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cut-off these functions are complicated. In the case of the proper-time cut-off 
[11] used here we have [9] 



r r A 2 

9 ( S ) = f^f*se* P {-s[e + *\} = ^E 2 



(7) 



and 



f(S,Q 2 ) = J— J dssj duexp{-s[k 2 + S 2 + u(l-u)Q 2 ]} 

Q 2 - 



16tt 2 



du Ei 



o 



+ u(l -u) 



(8) 



where Af is the quark cut-off, and the exponential integral is defined as 

oo _ . 

E n (x) = J dt^r. 

The one meson-loop gap equation (6) requires also the introduction of a reg- 
ulator for meson momenta. In other words, we have to regularize the di- 
vergent integral over d A Q. In Ref. [9] this was achieved by the substitution 

A b 

/ d 4 Q — > 7r 2 / dQ 2 Q 2 , where A b was the four- dimensional Euclidean meson 
o 

momentum cutoff. In the present study at finite temperatures, we employ the 
three-dimensional cutoff procedure, i.e., we make the replacement 



/ 



d 4 Q 



47T 



A 6 

duo J dqq 2 
o 



(9) 



where Q = (o>,q) and q — |q|. The form (9) is convenient for the implemen- 
tation of the boundary conditions satisfied by temperature Green's functions. 
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3 Finite temperature 

For calculations at finite temperature T we shall adopt the imaginary time 
formalism [12]. This can be done by making the following replacement in the 
quark momentum integrals 

Here F(k) = F(E, k) is an arbitrary integrand, and the sum runs over the 
fermionic Matsubara frequencies Ej = (2j + l)nT. The integral over the meson 
four-momenta should be also replaced by the sum of the form (10). In this case, 
however, the sum runs over the bosonic Matsubara frequencies u n = 2irnT. 
With this prescription we can turn to the calculation of the functions which 
are the finite temperature analogs of g(S) and f(S, Q 2 ). We find 




Fig. 1. Diagramatic representation of Eq. (6). The cross represents the first term, 
the one-quark-loop contribution corresponds to the second term, and meson-loop 
terms represent subsequent terms. The solid lines represent the quark propagator 
l/(D<D), the dased lines correspond to the meson propagators K of Eqs. (5), the 
external dased line represents scalar-isoscalar coupling, and the vertices follow from 
the form of (D^D). 
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T 



i 



— / du El 



A/ A / A^ 



Analogously, the inverse meson propagators become 



(12) 



K a (S,n,q,T) = (4NJ(S,n,q,T)(uj 2 n + q 2 + 4S 2 ) + a 2 m/S) * , 
K n (S,n,q,T) = (4NJ(S,n,q,T)(u J 2 n + q 2 ) + a 2 m/Sy 1 . (13) 

Finally, we can write the finite-temperature analog of Eq. (6): 



a 2 {S - m) - 8N c Sg(S, T) + ]T / g 2 x 



2/(5, 0, 0, T) + — n, g, T) (cu 2 + q 2 + 4S 2 )) 



+ 3 



dS 2 

2/(5, 0, 0, T) + n, q, T)(cu 2 + q 2 ) 



K a (S,n,q,T) 



K w (S,n,q,T)\=0. 



(14) 



If chiral symmetry is broken, then the above equation has a nontrivial solution 
for S. The quark condensate and S are related by the formula 



(qq) = - a 2 (S - m) , 



(15) 



which follows immediately from the fact that (qq) = 5T(&)/5m and Eq. (3). 



4 Low-temperature expansion in the chiral limit 



Before presenting our numerical results for (qq)r let us consider the low- 
temperature expansion. As shown by Gasser and Leutwyler [2], in the chiral 
limit the low-temperature expansion of the quark condensate has the form 

(^)T = (gg)o(i-^-^ + ...). (16) 



First, let us do the N c counting in this formula. Since F w ~ 0(yfN~^), subse- 
quent terms in the expansion are suppressed by 1/N C . Since our one-meson- 
loop calculation accounts for first subleading effects in the 1/N C expansion, we 
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can hope for reproducing only the T 2 term in Eq. (16). Further terms would 
require more loops. 



Using standard techniques [12], the sum over the bosonic Matsubara frequen- 
cies in Eq. (14) can be converted to a contour integral in the complex energy 
plane. By deforming this contour we collect all contributions from the singu- 
larities of the integrand, weighted with the thermal Bose distribution. At low 
temperatures, the dominant contribution comes from the lowest lying pion 
pole, and other singularities are negligible. Thus, the third term in (14) be- 
comes 



Af, 

3 JvI dqi 



> 1_ _3_ 

u 2 + q 2 2-7T 2 



dq q 



1 + 



(17) 



Writing Eq. (17) we have approximated the function f(M,n,q,T), appearing 
in the pion propagator, by its value at n = q = 0. For sufficiently large cutoff 
Aft, the integral over the thermal distribution function in (17) can be expressed 
by the Riemann zeta function ((2) = n 2 /Q. Thus, the final result for (17) is 
3A 2 /47r 2 + T 2 /2. Inserting the above result into the gap equation (14) we find, 
with m = 0, the following equality: 

SA 2 

h ( S , T)=a 2 - 8N c g(S, T) + — f + \T 2 = 0. (18) 



Eq. (18) defines implicitly the function S(T), which satisfies the equation 



dS dh(S,T)/dT 2 
dT 2 ~ ~ dh{S,T)/dS 



1QN, 



dg(S,T) 
OS 



(19) 



Here we have neglected the term dg(S,T)/dT 2 , since it is exponentially sup- 
pressed by the factor exp(— S/T). Furthermore, the leading-A c term on the 
right hand side of (19) can be rewritten using the relations [9] dg(S,T)/dS = 
-2Sf(S,0) and AN c f(S, 0) = F^/S 2 , where is the leading-^, piece of 
the pion decay constant. Collecting these equalities we arrive at dS/dT 2 = 
—S/(SF n ), which finally gives 



S(T) = S(0) 



1 - 



T 2 
8^ 



(20) 



Proportionality (15) implies that the above expression coincides (in the large 
7V C limit) with Eq. (16). Hence our method is consistent with a basic require- 
ment of chiral symmetry at the one-meson-loop level. 
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5 Results 



In the exact chiral limit the model has 3 parameters: a, Af, and A 6 . In this 
paper we fix arbitrarily A&/A/ = \. The remaining 2 parameters are fixed by 
reproducing the physical value of F n = 93MeV and a chosen value for (qq)o- 
For the case of m ^ we have an extra parameter, m, which is fitted by 
requiring that the pion has its physical mass. We compare results with meson 
loops to results with the quark loop only (A 6 = 0). Parameters for the two 
calculations are adjusted in such a way, that the values of F n , (qq)o, and m T 
are the same. 

The calculation of F n with meson loops, although straightforward, is rather 
tedious, so we do not present it here. The method has been presented in 
detail in Ref. [9,13]. The only difference in our calculation is that the three- 
momentum cut-off (9) rather than the four-momentum cut-off of Ref. [9] is 
used. 

Figure 2 shows the dependence of (qq) on T 2 . The solid line represents the 
case with meson loops. We note that at low temperatures the curve has a 
finite slope, as requested by Eq. (20). The slope is close to the leading-order 
Gasser-Leutwyler result (dotted curve). As explained earlier, the slopes would 
overlap in the large-A^ limit. This behavior is radically different from the 
case with quark loops only (dashed curve). In this case at low temperatures 
(qq)t — (<?<?) o ~ e~ M l T ', where M is the mass of the constituent quark. All 
derivatives of this function vanish at T = 0, and (qq) is flat at the origin. We 
can also see from the figure that the fall-off of the condensate is faster when 
the meson loops are included. In fact, for the parameters of Fig. 2 we have an 
interesting phenomenon. At T = 162MeV the condensate abruptly jumps to 0. 
There is a first-order phase transition, with a latent heat necessary to melt the 
quark condensate. Such a behavior is not present in the case of calculations 
without meson loops [4-6]. We note that with meson loops present the chiral 
restoration temperature is 162MeV, i.e. about 10% less than 176MeV of the 
quark-loop-only case. 

Figure 3 shows the same study, but for the physical value of m T . We note that 
now (qq) (solid line) is also flat at the origin, since the pion is no more massless, 
and at low T we have (qq)r ~ (qq)o ~ e~ m7r ^ T . Nevertheless, the region of this 
flatness is small, and at intermediate temperatures the curve remains close to 
the Gasser-Leutwyler expansion. We note again that meson loops considerably 
speed up the melting of the condensate compared to the case of quark loops 
only. However, there is no first-order phase transition such as in Fig. 2. Instead, 
we observe a smooth cross-over typical for the case of m ^ 0. 

The faster change of the quark condensate in our study is not surprising. It 
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is caused by the presence of light pions which are known to play a dominant 
role at low-temperatures [7]. The behavior of (qq) reflects this general feature. 
Concluding, we stress that the inclusion of meson loops in the NJL model 
qualitatively and quantitatively changes the results in comparison to the cal- 
culations at the quark-loop level. In particular, we find finite slope of (qq) vs. 
T 2 at the origin in the chiral limit, faster melting of the condensate, and lower 
chiral restoration temperature. 
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